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Abstract 

To  characterise  explosives  materials,  investigations  for  four  well  experimentally  characterised  energetic 
materials,  Cyclotrimethylenetrinitramine  (RDX),  Octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine  (HMX), 
Pentaerythritol  tetranitrate  (PETN)  and  Triamino-trinitrobenzene  (TATB)  were  conducted  using  quantum 
mechanics  and  analytical  methods.  Furthermore,  using  the  pressure  p  and  the  ratio  of  specific  densities, 
v/Vfl,  p-v  Hugoniot  diagrams  of  the  four  explosives  were  obtained.  Using,  the  Stine  equation  [Stine,  J. 
Energ.  Mater.,  1990]  and  more  recent  least  squares  fitting  algorithms,  new  constants  were  obtained,  to 
determine  the  detonation  velocity,  D ,  which  was  then  used  to  define  the  Raleigh  line.  The  ratio  of  specific 
heats,  y,  a  value  anywhere  between  2  and  3  and  usually  approximated  from  experimental  data,  was 
obtained  by  making  the  explosive  products  Hugoniot  tangent  to  the  Rayleigh  line  through  varying  y.  The 
effect,  sensitivity  and  importance  of  having  the  correct  value  of  y  were  demonstrated.  The  predicted 
parameters  namely  the  shock,  particle  and  detonation  velocities,  the  CJ  pressure  and  density,  the  ratio  of 
specific  heats  and  the  Hugoniot  diagrams  all  compared  quite  well  with  published  experimental  data. 

1.  Introduction 

To  better  design  energetic  materials  and  predict  their  energy  output  and  response  it  is  necessary  to 
characterise  these  materials  in  terms  of  their  energetics.  Indeed,  this  characterisation  demands  a  series  of 
experiments  ranging  from  slow  cook  off  of  the  material  to  high  velocity  flyer  plate  impact  tests  are 
necessary  to  obtain  parameters  such  as  the  detonation  velocity,  densities,  shock  and  particle  velocities  and 
critical  pressures.  These  experimental  parameters  are  required  to  obtain  the  Chapman- Jouguet  parameters 
that  define  the  point  at  which  the  detonation  occurs.  The  mean  objective  of  this  work  is  to  conduct 
numerical  experiments  to  examine  the  parameters  that  define  detonation  process  and  establish  a 
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completely  numerical  approach  by  using  quantum  mechanics  and  analytical  methods.  To  this  end,  four 
well  known  and  experimentally  characterised  energetic  materials  were  considered  and  the  results 
compared  with  experimental  values.  The  four  materials  chosen  are  RDX,  HMX,  PETN  and  TATB.  The 
detonation  velocity  is  one  of  the  key  performance  characteristics  of  energetic  materials  which  can  be 
determined  by  some  computer  codes  f  1-4].  Also,  pressure-specific  volume  relationship  called  Hugoniot 
curve  is  a  valuable  tool  to  understand  the  material  response  and  evaluate  their  performance.  A  treatment  of 
the  Hugoniot  curve  is  found  in  the  book  Explosives  Engineering  by  Cooper  [5].  Each  material  has  a 
unique  Hugoniot  curve  that  is  usually  obtained  from  a  series  of  experimental  tests  that  could  be  quite 
expensive  and  time  consuming.  In  this  work,  the  main  focus  was  on  the  characterisation  of  the  detonation 
of  the  energetic  materials  (including  the  Hugoniot)  using  the  numerical  simulations  technique.  This 
analysis  method  provided  a  valuable  way  to  understand  the  detonation  process  and  obtain  theirs 
corresponding  parameters.  The  computation  of  detonation  parameters  by  computer  codes  usually  requires 
the  heat  of  formation  AHf  and  the  density  of  energetic  materials;  therefore,  the  knowledge  of  heat  of 
formation  of  materials  remains  the  key  issue.  In  this  work,  the  hybrid  density-functional  theory  (DFT)  [6] 
approach  was  used  to  predict  heats  of  formation  in  the  solid  phase  and  detonation  properties  using 
quantum  chemistry  calculations. 

2.  Computational  methods  and  details 

All  calculations  on  structural  properties  of  compounds,  were  carried  out  in  the  framework  of  density- 
functional  theory  (DFT)  [7].  In  the  present  work,  B3LYP  hybrid  functional  [8H0]  was  used  in  conjunction 
with  6-31 1++G(3d,3p)  basis  sets.  The  geometry  of  molecules  were  optimized  using  gradient  techniques  at 
the  HF/STO-3G,  HF/6-31G,  B3LYP/6-311G,  and  B3LYP/6-311++G(3d,3p)  levels  without  symmetry 
constraints.  All  calculations  were  performed  using  the  Gaussian  03  rill  packages.  The  B3LYP  functional 
is  a  reliable  method  providing  accurate  results  when  compared  to  experimental  data  and  ab  initio  methods 

Q2]. 

In  this  study,  heats  of  formation  AfH^olid  were  determined  using  Hess’  law  of  constant  heat  summation 
[13]  given  by: 

HSub 


(1) 
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where  AfH^as  is  the  heat  of  formation  at  the  gas  phase  predicted  by  using  isodesmic  reaction  approach 

[4].  In  this  approach,  the  number  of  electron  pairs  and  chemical  bond  types  are  assumed  to  be  conserved 
in  the  chemical  reaction. 

The  heat  of  sublimation,  A  fHSub  ,  is  estimated  using  Politzer  approach  [14].  Politzer  et  al.  predict  heats  of 

sublimation  of  a  compound  by  considering  a  statistical  approach  where  various  parameters  are  mapped 
onto  isodensity  surfaces  of  the  isolated  molecule,  and  propose  an  empirical  expression  which  provides  the 
heat  of  sublimation  of  a  particular  compound  given  by: 

Ml  Sub  =  M2  +  A  (voi  )° 5  +  A  (2) 

where  A  is  the  molecular  surface  area.  It  is  the  collection  of  points  where  97%  of  the  electron  density  of 
the  isolated  molecule  is  present.  The  parameter  v  defines  the  balance  value  for  the  positive  and  the 
negative  potential  quantities  on  the  molecular  surface  A.  crfot  =  <r++<ri  represents  the  variance  of 

electrostatic  potential  values  on  the  molecular  surface  and  (3i_3  represent  the  fit  parameters  computed  from 
least-square  fitting  techniques.  The  method  proposed  by  Stine  [15]  to  obtain  the  detonation  velocity,  D, 
was  used.  A  least  squares  fitting  was  applied  to  minimise  the  sum  of  the  squares  of  the  residuals  formed 
by  the  difference  between  the  calculated  and  the  experimental  values  and  the  newer  fit  constants  were 
estimated.  The  Rayleigh  line  and  Hugoniot  curve  of  gas  products  were  obtained  using  the  simplest  theory 
discussed  by  Fickett  et  al.  [J_6].  The  model  assumes  that  the  products  of  reactions  are  instantaneously 
produced  by  the  shock.  The  conservation  of  mass,  momentum  and  energy  equations  were  used  and  the  gas 
products  of  reaction  were  assumed  as  an  ideal  gas.  The  four  materials  to  be  characterised  are  RDX,  HMX, 
PETN  and  TATB  and  their  chemical  structures  are  illustrated  in  Figure  1. 


RDX  HMX  PETN  TATB 


Figure  1.  Chemical  structure  of  molecules  in  the  studied  system. 
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3.  Results  and  discussion 

3. 1.  Gas -phase  heat  of formation 

In  this  work,  only  PETN  molecule  case  is  provided  as  an  example  for  gas-phase  heat  of  formation 
calculations.  The  gas-phase  heat  of  formation  of  PETN  was  estimated  using  isodesmic  reactions  given  by 

equation  3.  Using  the  experimental  NfHGas  values  of  Nitroglycerin  (NG)  and  Ethyl  nitrate,  A  fHGas  of 

PETN  was  calculated  as  follows: 


GTN  Ethyl  nitrate  PETN 


The  heat  of  reaction  of  this  equation  is  given  by: 

A# „„  =  ( )  -  (ArH Is, no 


I  A  ZJ® 

f  Gas, Ethyl  nitrate 


(3) 


(4) 


Using  the  experimental  heat  of  formation  of  the  molecules  making  up  the  PETN  and  available  in  NIST 
website1,  the  PETN  heat  of  formation  can  be  calculated  from  equation  4: 

Gas, PETN  ~  rxn  +  Gas, NG  +  ^  f  ^ Gas ,  Ethyl  nitrate  ~  ^f^Gas,H2  (5) 

Heats  of  formation  of  the  molecules  in  equation  3  are  listed  in  Table  1. 

Table  1.  Experimental  and  calculated  gas-phase  heats  of formation  of  the  molecules  making  up  the  PETN 


Molecules 

S/Hqos  ’  kcal/mol 

NG 

-67 

Ethyl  nitrate 

-37 

PETN 

-86 

h2 

0 

3.2.  Heats  of  sublimation 


1  www.  webbook.nist.  gov/chemistry 
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For  heats  of  sublimation,  only  PETN  case  is  studied  in  this  work.  Table  2  lists  some  similar  molecules  to 
PETN  which  have  been  used  to  determine  the  fit  parameters  Pi_3  shown  in  equation  2.  In  the  framework  of 
the  Politzer  approach,  and  the  heats  of  sublimation  A/Z^were  estimated.  The  experimental  values  of 

A HSub ,  of  the  molecules  used  in  the  fit  were  taken  from  the  website  www.  webbook.nist.  gov /chemistry. 
Politzer  parameters  shown  in  equation  2  are  A,  v  and  atotai  and  given  in  Table  2.  Figure  2  shows  the  plot  of 
predicted  A HSub  of  the  selected  molecules  similar  to  PETN.  The  curve  shows  a  good  correlation  between 
the  predicted  and  the  measured  data.  The  fitting  procedure  leads  to  Pi  =  2  x  10'4  kcal/mol.A4,  p2  =  2.104, 
and  p3  =  0.2374  kcal/mol.  Pi_3  in  equation  (2)  leads  to,  A HSub=  25.0kcal/mol  of  PETN.  Now  using 

equation  (1),  the  heat  of  formation  in  the  solid  phase  of  PETN  is  estimated  to  be  -  111  kcal/mol.  The 
experimental  value  as  reported  in  reference  [F7]  is  -128  kcal/mol,  the  relative  error  between  these  two 
values  is  about  15%.  This  deviation  can  be  possibly  reduced  by  using  extended  basis  set  in  the  hybrid 
functional  B3LYP. 


Figure  2.  Predicted  heats  of  sublimation  of  the  selected  molecules  versus  the  experimental  values 


Table  2.  Molecular  Properties. 


Names 

Chemical  structures 

A  (A2)  V 

^total 

(kcal/mol) 

A Hsub 

(kcal/mol) 

1, 1,1,3  -Tetranitro- 
2-methylpropane 

\  V 

'r\ 

o  ) - 

% 

- N 

y 

X° 

235.350  0.125 

10.963 

91.200 
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1,1,1,4-Tetranitr 

obutane 


2, 2,3,3- 

Tetranitrobutane 


Urea 


95.830  0.2472 


1, 1,1, 2,2- 

Pentanitropropane 


Ethyl  carbamate 


Trinitromethane 


N-Nitrosobis- 

(2,2,2- 

trinitroethyl)amine 


228.030  0.0596 


135.060  0.2167 


147.910  0.0378 


284.250  0.0521 


14.962 


12.394 


19.361 


12.561 


16.026 


16.571 


18.059 


99.600 


78.200 


98.600 


77.400 


76.300 


46.700 


23.326 
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N-Nitrobis-(2,2,2- 

trinitroethyl)amine 


.  U 
// 


OttttN 


295.280  0.0581 


13.891 


24.720 


3.3.  Determination  of  the  detonation  velocity 

Starting  with  the  empirical  method  proposed  by  Stine  [15]  to  obtain  the  detonation  velocity,  D ,  equation 
(6),  for  a  given  compound: 

D  —  D0  +  p  [cxnc  +  c2nN  +  c3n0  +  c4nH  +  c5Aj-HSolid  /  M  (6) 

Where  D0  is  a  constant,  nh  i  =  C,  N,  O ,  H ,  are  the  specific  number  of  atoms  in  the  molecule,  and  c/?  i  =  1 
to  5,  are  the  corresponding  characteristic  constants  for  the  carbon  (C),  nitrogen  (N),  oxygen  (O)  and 
hydrogen  (H)  atoms  in  the  molecule,  and  the  heat  of  formation  (A fH°Solid),  respectively,  p  is  the  density  and 

M  is  the  molecular  weight  of  the  material.  Equation  (6)  is  expressed  in  terms  of  a  set  of  known  and  well 
characterised  explosives,  for  which  experimental  detonation  velocities,  heats  of  formation  and  densities 
are  available.  Using  a  least  square  fitting  approach  to  equation  (6),  the  constants  D0  and  ct  could  be 
obtained.  These  constants  can  be  used  to  predict  the  detonation  velocities  of  new  compounds  that  have 
similar  chemical  composition. 

3.4.  Determination  of  constants  D0  and  ct 

The  base  set  of  explosives  compounds  selected  to  obtain  the  constants  D0  and  ct  are  exactly  the  same  five 
known  and  well  characterised  explosives  chosen  by  Stine  in  his  original  work  [151.  These  molecules  are 
BTNEU,  HMX,  TNETB,  ABH  and  Exp  ID.  Table  3  shows  the  equivalent  chemical  names  along  with  all 
the  necessary  basic  parameters  such  as  the  density,  experimental  detonation  velocity,  Dexp  and  molecular 
weight  that  are  required  in  equation  (6). 

Table  3.  Chemical  characteristics  of  BTNEU,  HMX,  TNETB,  ABH,  and  Exp  ID. 


Chemical 

Formula 

Compound 

Enthalpy  of 
Formation 
(kcal/mol) 

Density 

(g/cm3) 

Mol  Weight 
(g/mol) 

Dexp  km/s 

C24H6N14O24 

ABH 

116.3 

1.78 

874 

7.60 

c6h0n6o6 

BTF 

144.5 

1.86 

252 

8.49 

c5h6n8o13 

BTNEU 

-72.9 

1.86 

386 

9.01 
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c4h8n8o8 

HMX 

17.9 

1.89 

296 

9.12 

c6h6n6o14 

TNETP 

-118.5 

1.78 

386 

8.46 

c6h6n4o7 

Exp  ID 

-94.0 

1.55 

246 

6.85 

In  order  to  compute  the  detonation  velocity  D  of  RDX,  HMX,  PETN  and  TATB,  parameters  given  in 
Table  3  are  used  in  equation  (6)  and  then  a  least  squares  fitting  is  applied  to  minimise  the  sum  of  the 
squares  of  the  residuals  formed  by  the  difference  between  the  calculated  value  of  D  and  the  experimental 
value  Dexp,  and  the  D0 ,  ch  c2,  c3 ,  c4  and  c5  constants  are  estimated.  The  constants  D0  and  ct  (i  =  1  to  5) 
obtained  from  the  least  squares  fitting  of  equation  (6)  are  given  in  Table  4. 


Table  4.  New  data  set  of  D0  and  Cf  parameters  to  be  used  in  equation  (6). 


D0  (km/s) 

Cl 

C4 

<?5 

New  data 

2.96 

-7.039 

51.586 

69.034 

2.214 

0.463 

Stine 

3.69 

-13.85 

37.74 

68.11 

3.95 

0.6917 

The  least  square  fitting  produced  new  data  of  the  D0  and  ct  parameters  that  differ  from  Stine’s  coefficients. 
This  deviation  is  simply  due  to  the  fact  that  in  this  work  more  optimized  least  square  fitting  algorithms 
have  been  used. 

Table  5.  Detonation  velocities  using  Stine  and  new  constants  in  equation  (6). 


Compound 

DgXp  km/s 

D /Stine  km/s 

D  /New  data 

km/s 

%  of  error, 
Stine  data2 

%  of  error, 
New  data3 

ABH 

7.60 

7.63 

7.60 

0.3990 

0.0007 

BTF 

8.49 

8.50 

8.49 

0.1415 

0.0003 

BTNEU 

9.01 

8.95 

9.01 

0.6775 

0.0004 

HMX 

9.12 

9.02 

9.12 

1.0090 

0.0001 

TNETB 

8.46 

8.48 

8.46 

0.2300 

0.0006 

Exp  ID 

6.85 

6.86 

6.85 

0.1644 

0.0000 

From  Table  5,  the  calculated  %  error  using  new  data  is  lower  than  this  using  Stine’s  values.  For  this  raison 
the  named  constants  new  data  are  used  for  all  estimations  of  detonation  velocities.  The  value  of  each  term 
in  equation  (6)  has  been  determined,  now  it  is  possible  to  calculate  the  detonation  velocities  for  the 
selected  explosives  for  this  study  -  PETN,  RDX,  HMX  and  TATB.  The  experimental  heats  of  formation 
at  the  solid  states,  extracted  from  NIST  website4,  for  RDX,  HMX  and  TATB,  have  been  used  in  equation 
(6)  for  the  calculations  of  the  detonation  velocities  D.  Table  6  shows  a  comparison  between  the  calculated 


2  %  of  error  =  100*  \(Dexp-D/Stine)\/Dexp 

3  %  of  error  =  100*  |  (Dexp-D/New  Data)  \/Dexp ,  in  the  case  of  the  new  constants 

4  www. webbook.nist. gov/chemistry 
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and  measured  detonation  velocities  using  equation  (6)  and  very  good  agreement  is  obtained  between 
experimental  and  calculated  values  of  D.  Note  that  HMX  was  listed  as  a  molecule  in  the  Stine’s  set. 

Table  6.  Detonation  velocities  for  PETN,  RDX,  HMX  and  TATB 


Chemical 

Formula 

Compound 

Dexp,  km/s[-] 

D,  km/s 

c3h6n6o6 

RDX 

8.70 

8.73 

c4h8n8o8 

HMX 

9.12 

9.12 

c5h8n4o12 

PETN 

8.26 

8.29 

c6h6n6o6 

TATB 

7.76 

7.91 

3.5.  Hugoniot  curve  and  Rayleigh  line 

The  Rayleigh  line  and  Hugoniot  curves  of  gas  products  are  obtained  using  the  simplest  theory  and  are 
discussed  in  many  shock  physics  texts  such  as  Fickett  et  al.  [16].  The  model  assumes  that  the  products  of 
the  reaction  are  instantaneously  produced  by  the  shock.  Using  the  conservation  of  mass,  momentum  and 
energy  equations,  and  the  equation  for  Rayleigh  line  is: 

■^-0  (7) 

K-V) 

Which  is  expressed  in  terms  of  the  specific  volume  v  =  Up ,  pressure  and  detonation  velocity 

For  instance,  an  ideal  gas  is  characterised  with  a  constant  heat  capacity  and  a  constant  heat  for  its 

complete  reaction  q  =  AH°r  and  the  Hugoniot  curve  that  is  expressed  as, 


f 


P  2 

—  +  JU 

Po 


- P 

\vo 


-1  +  JU4  —  fJL 


2  Tkq  _ 


Povo 


=  0 , 


(8) 


where  ju  =  (y  -  1)1 (y  +  1)  and  y=  Cp/Cv  is  the  ratio  of  specific  heats.  X  is  the  degree  of  the  chemical 
reaction,  changing  from  0  for  no  reaction  to  1  for  a  complete  reaction.  In  this  work,  it  is  assumed  that  the 
reaction  is  complete,  then  X=\. 


This  is  the  equation  of  a  rectangular  hyperbola  in  the  p/p0  -  v/v0  plane,  centered  at  the  point  v/v0  =  p  and 

p/po  =  -  fC. 


If  the  Hugoniot  curve  equation  (8)  is  solved  simultaneously  with  the  Rayleigh  line  equation  (7),  the 
intersection  of  these  two  equations  defines  the  CJ  point  at  the  detonation  process,  q  represents  the  heat  of 
detonation;  it  can  be  estimated  using  quantum  mechanical  method  [J_8].  In  this  work,  the  heat  of 
detonation  q  and  the  density  p0  values  are  extracted  from  references  [3,  19,  20]  and  used  in  equations  (7) 
and  (8)  to  plot  Hugoniot  curves  of  gas  products  and  Rayleigh  lines.  The  line  tangent  at  the  CJ  point  (pCJ, 
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Vcj)  is  the  Rayleigh  line.  The  slope  of  line  tangent  at  the  CJ  point  and  the  Rayleigh  line  are  equal.  The 
analytical  results  for  the  four  explosives,  at  various  corresponding  values  of  y,  are  presented  in  Figure  3. 

3. 6.  The  Gamma,  y  effect 

The  detonation  process  of  an  explosive  the  reaction  occurs  immediately  and  this  implies  that  once  the 
unreacted  material  passes  through  the  detonation  wave.  It  is  assumed  that  it  is  instantaneously  converted 
to  a  new  material  and  as  discussed  previously  the  final  state  is  determined  by  the  intersection  of  equation 
(8)  and  the  Raleigh  line  equation  (7).  The  correct  value  of  y  which  is  the  ratio  of  the  specific  heat  at 
constant  pressure  to  that  at  constant  volume  ( CPICV )  is  important  because  it  influences  significantly  the 
tangent  of  the  Hugoniot  energy  and  as  a  result  will  influence  the  CJ  point.  As  pointed  out  by  Cooper  [5], 
an  explicit  value  of  y  is  not  available  but  experiments  have  shown  that  it  ranges  between  2  and  3.  Many 
texts  provide  approximate  values  for  y.  However,  in  this  work  the  idea  was  also  to  explore  the  influence  of 
y  on  equation  (8).  For  given  values  of  p  and  v,  y  was  varied  between  2  and  3  for  the  four  explosives. 
Figure  3  shows  the  effect  of  y  on  equation  (8)  and  its  relation  to  intersecting  the  Raleigh  line,  for  the  four 
molecules.  In  the  case  of  PETN,  when  y  is  equal  to  2.2,  the  equation  (8)  intersects  the  Raleigh  line  at  two 
points  A  and  B  and  as  y  is  gradually  increased  (around  2.55)  at  some  point  it  is  tangent  to  the  Raleigh  line 
and  if  further  increased  (3.0)  it  stays  above  the  Raleigh  line.  Similar  analysis  is  observed  for  HMX,  RDX 
and  TATB.  In  this  way  the  CJ  point  for  the  four  explosives  were  obtained  and  these  are  given  in  Table  7 
and  the  results  reported  in  Figure  3.  The  values  of  y  compare  very  well  with  published  experimental  data. 
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Table  7.  Experimental  and  calculated  properties  of  explosives. 


r\A  ™ 

r 

Pcj  (GPa) 

pu  pxg/m  ) 

Theo 

Exp  L-J 

Theo 

ExpL1J 

RDX 

1800 

2.83 

2.98 

36.45 

34.70 

HMX 

1900 

2.89 

3.00 

40.60 

39.30 

PETN 

1770 

2.55 

2.64 

34.50 

33.50 

TATB 

1895 

3.00 

2.72 

29.23 

31.50 

Figure  3.  Hugoniot  curves  to  demonstrate  the  y  effect  relative  to  the  Raleigh  lines  for  the  four  explosives. 

4.  Conclusions 

Through  the  investigation  conducted  in  this  for  four  well  experimentally  characterised  energetic  materials, 
RDX,  HMX,  PETN  and  TATB  using  quantum  mechanics  calculations  coupled  with  analytical  methods,  it 
is  now  possible  to  estimate  some  parameters  such  as  the  detonation  velocity  and  the  Chapman- Jouguet 
pressure  and  density  that  define  the  detonation  of  a  CHNO  based  explosives.  .  The  detonation  velocities  of 
explosives  were  obtained  using  Stine’s  equation.  However,  new  and  seemingly  more  accurate  constants 
were  obtained  for  this  equation  using  instead  a  more  optimized  algorithm  for  a  least  squares  fitting.  With 
these  new  accurate  constants  leading  to  precise  predicted  detonation  velocities  when  they  were  compared 
to  experimental  values.  The  specific  heat  ratio,  a  parameter  found  in  the  Hugoniot  energy  equation  for  the 
products,  and  the  Chapman- Jouguet  pressures  and  specific  densities  were  obtained  for  the  studied 
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explosives  using  the  point  where  the  Raleigh  line  tangents  to  this  Hugoniot  diagram.  The  values  obtained 
for  Chapman- Jouguet  parameters  and  the  specific  heat  ratios  compared  very  well  with  published 
experimental  values  for  all  four  explosives.  A  sensitivity  analysis  of  the  specific  heat  ratio  was  examined 
and  it  was  shown  that  the  tangent  that  energy  Hugoniot  makes  with  the  Raleigh  line  is  quite  sensitive  and 
care  needs  to  be  taken  in  its  determination  since  the  Chapman- Jouguet  parameters  are  determined  using 
this  tangent  point.  To  conclude,  the  potential  use  of  the  molecular  modelling  approach  together  with  the 
analytical  methods  shows  that  now  a  CHNO  based  explosives  can  be  numerically  characterised  and 
parameters  such  as  CJ,  specific  heat  ratio,  y,  pCj  and  detonation  velocities  can  be  obtained  quite  accurately 
prior  to  experiments.  This  approach  should  be  also  as  supporting  information  to  suggest  a  minimum 
number  of  or  just  confirmation  experiments  if  necessary. 
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